An unfitted discontinuous Galerkin scheme for conservation laws on evolving surfaces

Motivated by considering partial differential equations arising from conservation laws posed on evolving surfaces, a new numerical method for an advection problem is developed and simple numerical tests are performed. The method is based on an unfitted discontinuous Galerkin approach where the surface is not explicitly tracked by the mesh which means the method is extremely flexible with respect to geometry. Furthermore, the discontinuous Galerkin approach is well-suited to capture the advection driven by the evolution of the surface without the need for a space-time formulation, back-tracking trajectories or streamline diffusion. The method is illustrated by a one-dimensional example and numerical results are presented that show good convergence properties for a simple test problem

The interface for functions in the dune-functions module

The dune-functions dune module introduces a new programmer interface for discrete and non-discrete functions. Unlike the previous interfaces considered in the existing dune modules, it is based on overloading operator(), and returning values by-value. This makes user code much more readable, and allows the incorporation of newer C++ features such as lambda expressions. Run-time polymorphism is implemented not by inheritance, but by type erasure, generalizing the ideas of the std::function class from the C++11 standard library. We describe the new interface, show its possibilities, and measure the performance impact of type erasure and return-by-value

Function space bases in the dune-functions module

The dune-functions Dune module provides interfaces for functions and function space bases. It forms one abstraction level above grids, shape functions, and linear algebra, and provides infrastructure for full discretization frameworks like dune-pdelab and dune-fem. This document describes the function space bases provided by dune-functions. These are based on an abstract description of bases for product spaces as trees of simpler bases. From this description, many different numberings of degrees of freedom by multi-indices can be derived in a natural way. We describe the abstract concepts, document the programmer interface, and give a complete example program that solves the stationary Stokes equation using Taylor-Hood elements

Co-assembly of chitosan and phospholipids into hybrid hydrogels

Novel hybrid hydrogels were formed by adding chitosan (Ch) to phospholipids (P) self-assembled particles in lactic acid. The effect of the phospholipid concentration on the hydrogel properties was investigated and was observed to affect the rate of hydrogel formation and viscoelastic properties. A lower concentration of phospholipids (0.5% wt/v) in the mixture, facilitates faster network formation as observed by Dynamic Light Scattering, with lower elastic modulus than the hydrogels formed with higher phospholipid content. The nano-porous structure of Ch/P hydrogels, with a diameter of 260±20 nm, as observed by cryo-scanning electron microscopy, facilitated the penetration of water and swelling. Cell studies revealed suitable biocompatibility of the Ch/P hydrogels that can be used within life sciences applications

The Dune framework: Basic concepts and recent developments

This paper presents the basic concepts and the module structure of the Distributed and Unified Numerics Environment and reflects on recent developments and general changes that happened since the release of the first Dune version in 2007 and the main papers describing that state Bastian etal. (2008a, 2008b). This discussion is accompanied with a description of various advanced features, such as coupling of domains and cut cells, grid modifications such as adaptation and moving domains, high order discretizations and node level performance, non-smooth multigrid methods, and multiscale methods. A brief discussion on current and future development directions of the framework concludes the paper

Cell-scale degradation of peritumoural extracellular matrix fibre network and its role within tissue-scale cancer invasion

Local cancer invasion of tissue is a complex, multiscale process which plays an essential role in tumour progression. Occurring over many different temporal and spatial scales, the first stage of invasion is the secretion of matrix degrading enzymes (MDEs) by the cancer cells that consequently degrade the surrounding extracellular matrix (ECM). This process is vital for creating space in which the cancer cells can progress and it is driven by the activities of specific matrix metalloproteinases (MMPs). In this paper, we consider the key role of two MMPs by developing further the novel two-part multiscale model introduced in [33] to better relate at micro-scale the two micro-scale activities that were considered there, namely, the micro-dynamics concerning the continuous rearrangement of the naturally oriented ECM fibres within the bulk of the tumour and MDEs proteolytic micro-dynamics that take place in an appropriate cell-scale neighbourhood of the tumour boundary. Focussing primarily on the activities of the membrane-tethered MT1-MMP and the soluble MMP-2 with the fibrous ECM phase, in this work we investigate the MT1-MMP/MMP-2 cascade and its overall effect on tumour progression. To that end, we will propose a new multiscale modelling framework by considering the degradation of the ECM fibres not only to take place at macro-scale in the bulk of the tumour but also explicitly in the micro-scale neighbourhood of the tumour interface as a consequence of the interactions with molecular fluxes of MDEs that exercise their spatial dynamics at the invasive edge of the tumour

Structured models of cell migration incorporating molecular binding processes

The dynamic interplay between collective cell movement and the various molecules involved in the accompanying cell signalling mechanisms plays a crucial role in many biological processes including normal tissue development and pathological scenarios such as wound healing and cancer. Information about the various structures embedded within these processes allows a detailed exploration of the binding of molecular species to cell-surface receptors within the evolving cell population. In this paper we establish a general spatio-temporal-structural framework that enables the description of molecular binding to cell membranes coupled with the cell population dynamics. We first provide a general theoretical description for this approach and then illustrate it with two examples arising from cancer invasion

Toward transient finite element simulation of thermal deformation of machine tools in real-time

Finite element models without simplifying assumptions can accurately describe the spatial and temporal distribution of heat in machine tools as well as the resulting deformation. In principle, this allows to correct for displacements of the Tool Centre Point and enables high precision manufacturing. However, the computational cost of FE models and restriction to generic algorithms in commercial tools like ANSYS prevents their operational use since simulations have to run faster than real-time. For the case where heat diffusion is slow compared to machine movement, we introduce a tailored implicit–explicit multi-rate time stepping method of higher order based on spectral deferred corrections. Using the open-source FEM library DUNE, we show that fully coupled simulations of the temperature field are possible in real-time for a machine consisting of a stock sliding up and down on rails attached to a stand